Computer simulations of the mechanism of thickness selection in polymer crystals 
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In this paper I describe the computer simulations that I have performed to criticaUy examine 
the Lauritzen-Hoffman and the Sadler-Gilmer theories of polymer crystallization. In particular, I 
have computed the free energy profile for nucleation of a new crystalline layer on the growth face 
to compare with that assumed by the Lauritzen-Hoffman theory, I have analysed the mechanism 
of thickness selection in a multi-pathway model in which some of the constraints in the Lauritzen- 
Hoffman theory are relaxed, and I have re-examined the model used by Sadler-Gilmer. These 
investigations have lead to a mechanism of thickness selection of lamellar polymer crystals that 
differs from the two theories that I set out to examine. 



I. INTRODUCTION 



In 1957 Andrew Keller reported that polyethylcr 



1 



formed chain-folded lamellar crystals from solution 
This discovery was followed by the confirmation of 
the generality of this morphology — lamellar crystals are 
fornied on crystallization from both solution and the 
meltfi for a wide variety of polymers — and the basic 
phenomenological laws desojibing such properties as the 
thickness and growth rate.EpQ In particular, the crystal 
thickness, I, has beenpfound to be inversely proportional 
to the supercoolingjQ'u which is interpreted as resulting 
from / being slightly larger than /min, the minimum thick- 
ness for which a lamellar crystal is stable with respect to 
the solution or melt, i.e. I = Imin + 51, where 51 is small. 

Surprisingly, however, no theoretical consensus has yet 
been reached as to the mechanism of this seemingly sim- 
ple behaviour. In particular, two of the most well-known 
theoriea-ppthe Lauritzen-Hoffman (LH) surface nucleation 
theorjatELand the Sadler-Gilmer (SG) entropic barrier 
mode]ll3il3 — present very different explanations of thick- 
ness selectionJij Of course, in such a situation, one would 
like to determine which of the theories, if any, is closest to 
the truth. There are two aspects to such a task. Firstly, 
the predictions of the theories should be critically com- 
pared with experimental results. In the case of polymer 
crystallization both the LH and SG theories are able to 
reproduce the basic behaviour: the observed temperature 
dependence of the thickness and the growth rate. Addi- 
tionally, Hoffman and coworkers have further developed 
the surface nucleation approach in order to explain some 
of the more detailed behaviour of crystallizing polymers, 
for example the regime transitions in the growth rate.0 
However, this comparison does not conclusively favour 
one of the theories. This situation illustrates the fact 
that although consistency with experiment is an impor- 
tant first hurdle for any theory, it does not automatically 
imply the correctness of a theory. There may be a number 
of different ways of generating a particular experimental 
law. Furthermore, the number of parameters in a com- 
plex theory may give the theory sufficient plasticity to fit 
a wide variety of scenarios. 

Secondly, it is important that the assumptions of a 



theory, particularly those about the microscopic mech- 
anisms, are critically examined. However, in the case 
of polymer crystallization this task is very difficult to 
achieve experimentally. By addressing this gap, com- 
puter simulations can potentially play an important role 
in this field. Such simulations could range from examin- 
ing simple models to performing realistic atomistic sim- 
ulations of the crystal growth process. The former could 
allow the effects of relaxing some of the theoretical as- 
sumptions to be determined and the latter could provide 
a detailed molecular picture of the growth process. In- 
deed, there has been an increasing Jiupber of computa- 
tional studies pursuing these aims.ll3~t2 In this paper I 
will review my efforts in this directionLj~Ej and hope to 
illustrate the positive role that computer simulations can 
play in helping to understand polymer crystallization. In 
particular, the aim of my simulations has been to criti- 
cally examine the LH and SG theories. 



II. FREE ENERGY PROFILES 

In the LH theory the growth of a new layer is modelled 
as the deposition of a succession of stems (straight sec- 
tions of the chain that traverse the growth face) along the 
growth face from an initial nucleus, where the length of 
each stem is the same as the thickness of the lamella. The 
inset of Figure |l| illustrates the geometry of this mech- 
anism. To analyse the kinetics of growth, a thermody- 
namic description of the nucleation and growth of a new 
layer is first required. The free energy of a configuration 
with A"stom complete stems is taken to be 

^(A^stcm) = 26/a + 2(iVstcm - l)ahaf - N.tc^ablAF, 

(1) 

where a and b are the width and depth of a stem, I is 
the thickness of the lamella, a is the lateral surface free 
energy, af is the fold surface free energy, and AF is the 
free energy of crystallization. The first term corresponds 
to the free energy of the two lateral surfaces created on 
the deposition of the first stem and is proportional to I. 
The second term is the free energy of the new fold sur- 
face created on the deposition of subsequent stems. It is 



then assumed that at the barrier between configurations 
with different numbers of stems all the new surfaces have 
been created and that a fraction V? of the free energy of 
crystallization is released. This then gives the LH free 
energy profile that is illustrated in Figure m 




FIG. 1. The free energy profile for the nucleation and 
growth of a new layer assumed by the LH theory. The inset 
is a schematic representation of a configuration with three 
stems deposited. 

From this free energy profile, S{1), the flux over the 
barrier, can be obtained. The observed crystal thickness 
is then taken to correspond to the average 
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lS{l)dl. 



(2) 



This average thickness is close to the value of I at the 
maximum in S{1), which in turn is close to, but slightly 
above ?min, thus reproducing the observed behaviour of 
I. The maximum in S{1) is the result of two compet- 
ing factors. The free energy barrier for deposition of the 
first stem increases with I, thus making the growth of 
thick crystals prohibitively slow. However, as Imin is ap- 
proached from above, the thermodynamic driving force 
for crystallization goes to zero. 

It is important to note that by integrating over I, Equa- 
tion (0) assumes that there are crystals with all values 
of I greater than ^min which all grow with constant thick- 
ness and contribute to the average I. Those crystals with 
a thickness close to the maximum in S{1) dominate this 
ensemble and contribute more to Equation ^ because 
of their rapid growth. As was realized by Frank and 
Tosi,L3 the results of experiments where the temperature 
is changed during crystallization argue against such an 
ensemble. The temperature jumps give rise to steps on 
the lamellae, showing that a crystal need not necessarily 
grow at constant thickness.ESEZl 

We will come back to this issue later, but in this sec- 
tion we focus on the LH free energy profile. In particular, 
we compare this theoretical profile with-pnes computed 
from simulations of a simple polymerB3 In our model 



the polymer is represented by a self-avoiding walk on a 
simple cubic lattice. There is an attractive energy, -e, 
between non-bonded polymer units on adjacent lattice 
sites and between polymer units and the surface, and an 
energetic penalty, eg, for kinks (or 'gauche bonds') in the 
chain. The parameter eg determines the stiffness of the 
chains. In our simulations we have included a surface 
which represents the growth face of a polymer crystal. 
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FIG. 2. Free energy profiles for the formation of a target 
crystal with 5 stems of length 40 units and adjacent reentry 
of the folds. In (a) the profiles have been calculated from 
simulation, the labels give the temperature, and example 
configurations along the pathway have been illustrated.. In 
(b) the profile has been calculated at T = 2.75 ek~ using 
Equation (fBl) for pathways which allow one (1) or two (2) 
incomplete stems, as labelled, eg = 4e. 

To follow the crystallization of the polymer on the sur- 
face, we need to define an order parameter which de- 
termines the degree of crystallinity. We use A'xtai, the 
largest fragment of the polymer with the structure of a 
target crystalline configuration. In our case, we examine 
the crystallization of a 200-unit chain into a structure 
with 5 stems of length 40 units. In order to compare with 
the theoretical profiles we have to constrain the other 
N — A'xtai units in the chain to be disordered. The simu- 



lations-were carried out using configurationalpbias Monte 
CarlOjEa and the umbrella sampling techniquecj was used 
to calculate the free energy profiles. 

The free energy profiles that we obtained are shown in 
Figure 0a. They show the expected temperature depen- 
dence: at low temperature the crystal is most stable and 
at high temperature the disordered state is most stable. 
Note that the value of A'xtai for the disordered state is 
non-zero, because the disordered polymer is adsorbed on 
the surface. The adsorbed polymer is bound to have some 
short straight sections that qualify as crystalline by the 
definition of iVxtai- The free energy profiles also have a 
sawtooth structure resembling that of the theoretical pro- 
file. The barriers occur immediately after the previous 
stem has been completed, and correspond to the forma- 
tion of a new fold. They are followed by a monotonic de- 
crease in energy as this new stem grows to completion. In 
the language of the LH theory ^(A'^Btcm -^ ^stcm -I- 1) « 
for iVstom ^ 2. However, there is no feature in the sim- 
ulation profiles that corresponds to the formation of the 
first fold. This is because the initial nucleus is not a sin- 
gle stem, but two stems connected by a fold that grow 
simultaneously. Such a possibility had previously been 
suggested by Point .E3 

Confirmation of a two-stem nucleus comes from a sim- 
ple model calculation of the free energy profile. We can 
write the free energy as 

A(7V,tai) = ^coii(Af - iV,tai) + fcT^ exp {-E^te.l/kT) , 

(3) 

where the sum is over all possible crystalline configura- 
tions which are A'xtai units long, ii^xtai is ttiG energy of 
the crystalline configuration, and j4coii is the free energy 
of an ideal two-dimensional coil. The resulting profile 
is very similar to the simulation profile (Figure pp). In 
particular, there is no feature due to the formation of 
the first fold. However, when we force the initial nu- 
cleus to be a single stem by restricting the sum in the 
above equation to only those crystalline configurations 
with one incomplete stem, a free energy barrier associ- 
ated with the formation of the first fold appears. The 
reason for the preference for a two-stem nucleus is sim- 
ply energetic. For TVxtai > 4eg/e-|-2 the two-stem nucleus 
is lower in energy because of the interaction between the 
two stems. Our simulations were performed on a surface 
that was infinite. Whether a two-stem nucleus would be 
expected, when, as with a lamellar crystal, the thickness 
of the growth face is finite, depends upon how this critical 
size compares to the thickness of the lamella. 

It can be clearly be seen from Figure Ijb that the two- 
stem nucleus significantly reduces the nucleation barrier. 
In particular, it will no longer be proportional to I. This 
has significant implications for the LH theory given the 
key role played by this initial free energy barrier in con- 
straining I to a value close to Imm- 

Before we move on we should make a number of com- 
ments. First, the polymer model is very simple, and al- 



though there is no obvious reason why the thermody- 
namic reasons behind the two-stem nucleus should not 
also apply to a real polymer, there may be factors that 
are not included in our model that come into play. 

Second, the profiles refiect our choice of order param- 
eter. As we monitor crystallization unit by unit, during 
the growth of the first two stems the lateral surface en- 
ergy is paid for at the same time as the free energy of 
crystallization is released. Therefore, in this size range 
^ is effectively equal to 1, albeit with the possibility of a 
two-stem nucleus. Hoffman, however, advocates a ^=0 
version of the LH theory — it has the advantage that it 
avoids a '(5Z-catastrophe' (a divergence of the lamellar 
thickness) at large supercooling — in which he postulates 
that prior to crystallization an aligned physisorbed state 
is formed that has lost its entropjy but not yet gained 
the free energy of crystallization.B Such a state cannot 
occur in our lattice model because there is no difference 
in the interaction with the surface for a disordered chain 
adsorbed on the surface and a crystalline layer. In an 
off-lattice model the interaction energy for the crystalline 
layer would be greater because the stems would fit into 
the grooves provided by the steins of the previous layer. 



FIG. 3. An example crystalline configuration during the 
growth of a new layer on the surface of the growth face in the 
multi-pathway model. The arrows indicate the five possible 
moves at the next step in the kinetic Monte Carlo simulation. 
The dotted lines show the edges of the growth face. 

Third, a good order parameter must pass continuously 
through intermediate values when the system goes be- 
tween two states. However, one can imagine a number 
of mechanisms by which this criterion for Axtai is bro- 
ken. For example, in a realistic simulation of the sur- 
face crystallization of a long alkane into a once-folded 
configuration, the chain first formed non- adjacent crys- 
talline stems connected by a loose fold which then came 
together by-, the propagation of a defect through one of 
the stems.ta Another possibility that has been observed 
in simulations is the formation of crystallites in different 
portions of a chaittptiiat subsequently coalesce to form a 
single crystallitctSta 



III. A MULTI-PATHWAY MODEL 



In the previous section, in order to compare the LH 
free energy profile with those from simulation, we had 
to constrain the N — N^tsd units not having the target 
structure to be disordered. If we had not done this, at 
temperatures where the crystal is most stable the rest of 
the chain would have formed a crystalline configuration 
with stem lengths different from the target configuration. 
This naturally raises questions about the LH assumption 
that the stems in a new layer must all have the same 
thickness as the previous layer. In this section, we ex- 
amine the effects of relaxing some of the LH assumptions 
by studying a model in which the stems growpuaijt by 
unit and the length of a stem is unconstrained. EirE3 We 
term it a multi-pathway model because it can take into 
account the many possible ways that a new crystalline 
layer can form. 

This idea is not new.r-Jrank and Tosi,E3 PriceEJ and 
Lauritzen and Passagliaa considered models jKihere the 
stem length is not alwajys constant, and Point ,t3 and Di- 
Marzio and GuttmanE^ studied models where the stems 
could grow unit by unit. All these studies were performed 
at a time when computational resources were much less, 
so approximations and simplifications had to be made in 
order to render the models tractable. The natural way 
to solve such problems, though, is through the use of 
computational techniques, such as kinetic Monte Carlo. 
However, the only applications of computationaLmethods 
to this problem were in a short note by PointE3 and the. 
continuation of this work in the PhD thesis of Dupire.t^l 
Some of the results presented in these earlier studies are 
similar to those we report here. 

In our model we grow a single new crystalline layer 
by the successive growth of stems across a surface that 
represents the growth face of a polymer crystal. The 
polymer interactions are the same as used in the previ- 
ous section, and we only model the crystalline portion 
of the polymer explicitly — the rest is assumed to behave 
like an ideal coil. An example configuration is illustrated 
in Figure along with possible changes of configuration. 
These changes can only occur at the ends of the crys- 
talline portion, and are selected using the kinetic Monte 
Carlo algorithm, in which a move is chosen with a prob- 
ability proportional to the rate for that process. 

First, we shall examine the effect of the initial nu- 
cleus on the thickness of the layers grown. If the stem 
lengths are unconstrained and the initial nucleus is a sin- 
gle stem, one might imagine that one way of reducing the 
large initial free energy barrier in Figure n (and achiev- 
ing faster initial growth) would be for the stem length 
to increase gradually to its average value as crystalliza- 
tion progresses. For this pathway, the lateral surface free 
energy is paid for 'in installments' rather than all ini- 
tially. This is exactly what we observe when we force 
the initial nucleus to be a single stem by only allowing 
growth from one end of the crystalline portion of the 



chain (Figure ^a). When a double-stem nucleus is al- 
lowed the initial growth is very different because there is 
now no longer a large initial free energy barrier to cir- 
cumvent. The most important thing to note from these 
results is that, contrary to the LH theory, the thickness 
of the inital nucleus does not determine the thickness of 
the layer. Further confirmation of this can be obtained 
when we examine the growth from initial seed crystals. 
Whatever the thickness of the initial seed the thickness 
of the growing crystal converges to the same value (Fig- 
ure Eb). This implies that the thickness of a crystalline 
layer must be determined by factors which are operating 
on the deposition of each stem and not those specific to 
the initial stems. 
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FIG. 4. The dependence of tlie average stem length on 
the distance of the stem from (a) the initial nucleation site, 
and (b) the centre of an initial crystal seed for growth on 
an infinite surface at T = 2.75 efc~^. In (a) growth starts 
with a single polymer unit on the surface and we consider 
the cases where growth is allowed at one end or both ends 
of the crystalline configuration. In (b) the crystal seeds are 
3 stems wide; the lengths of the stems in the seeds are as 
labelled, eg = 8e. 

To determine what these factors might be, in Figure 
we show how the thickness of a new layer depends on 



temperature. First, it is immediately obvious that the 
thickness of a new layer is not necessarily the same as 
that of the growth face. Second, all the curves increase 
as the temperature approaches T„i, the melting or dis- 
solution temperature, because of the rise of Imm- Third, 
the thickness also increases at low temperature, in this 
instance because it becomes increasingly difficult to scale 
the free energy barrier for forming a fold and so on aver- 
age the stems continue to grow for longer. However, this 
rise is checked by the thickness of the growth face. It is 
unfavourable for the polymer to overhang the edge of the 
growth face because these units do not interact with the 
surface. 
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FIG. 5. The temperature dependence of the average stem 
length in a new crystalline layer for growth of a single layer 
on growth faces of different thickness, as labelled. 

Figure Isl only describes the growth of a single layer. 
However, as the thickness of the new layer is not gen- 
erally the same as the thickness of the growth face, one 
needs to consider the addition of a succession of layers. 
If we assume that all the variations in the stem length 
within a layer are annealed out before a new layer begins 
to grow, this can be achieved using Figure ga, in which 
we have plotted for a single temperature the thickness of 
the new layer against the thickness of the growth face. By 
following the dotted lines one can see what would hap- 
pen for growth on a growth face that is 50 units thick: 
the first layer is 36 units thick, the second 28, the third 
23, . . . Thus, the thickness converges to the value I** at 
which the curve crosses y = x, i.e. to the point where the 
thickness of the new layer is the same as the previous, 
and then the crystal continues to grow at this thickness. 
The mapping represented in Figure Ha is a fixed-point 
attractor. 

A similar picture emerges if we explicitly perform sim- 
ulations of multi-layer growth. Figure [t] shows a cut 
through a typical configuration that results. Within 5-10 
layers the thickness of the crystal converges to its steady- 
state value I** and then growth continues at that thick- 
ness. 
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FIG. 6. The dependence of the average stem length in 
a new crystalline layer on the thickness of the growth face 
for the growth of a single layer at (a) T — 2.75 efc~^ and (b) 
T — 3.375 efc~^. The dotted lines show how the thickness 
changes on addition of successive layers to a 50-unit thick 
surface. The inset in (a) shows the growth rate for the new 
layer as a function of the thickness of the growth face. 



The mechanism of thickness selection that occurs in 
our multi-pathway is at odds with the LH theory. It 
shows that it is inappropriate to compare the growth 
rates of crystals of different thickness because the thick- 
ness has only one dynamically stable value for which 
growth at constant thickness occurs. The ensemble of 
crystals assumed by Equation (0) is fictitious. Further- 
more, the growth rate of a new layer slows down as /** is 
approached from above (inset of Figure ||a). However, 
we shottld^note that in some of the multiple-pathway 
studiescjH'EJ mentioned earlier, it was realized that sta- 
ble growth can only occur at the one thickness where a 
new layer has the same thickness as the previous. Since 
then this insight has for the most part been neglected. 

To analyse the reasons for the dynamical convergence 
of the thickness to the value V* we examine how the 
probability distributions for the stem length depend on 
the thickness of the growth face (Figure @). /min places 



one constraint on the stem length; only a small fraction 
of the stems can be shorter than ^min if the layer is to be 
thermodynamically stable. The thickness of the growth 
face places the second constraint on the stem length; it 
is energetically unfavourable for the polymer to extend 
beyond the edges of the growth face. There is also a 
third weaker kinetic constraint on the stem length. At 
every step there is always a certain probability that a 
fold will be formed. Therefore, even in the absence of 
the second constraint, i.e. an infinitely thick growth face, 
the probability distribution will decay exponentially to 
zero at large stem length (Figure ||a). Although, this 
effect preveiitSj4ljie thickness from ever diverging in a Sl- 
catastrophej2j'Eil it does not stop the thickness becoming 
very large. 

When the growth face is significantly thicker than ^min 
there is a range of stem lengths between Zmin and the 
thickness of the growth face that are viable, and there- 
fore the new layer will be thinner than the previous layer. 
However, as the thickness of the growth face decreases, 
the probability distributions of the stem length becomes 
increasingly narrow and the difference in probability be- 
tween the stem length being greater or less than the sur- 
face thickness diminishes. Finally, at I** , as the thickness 
of the growth face approaches /min, the probability dis- 
tribution become symmetrical about the surface thick- 
ness and the thickness of the new layer becomes equal to 
the thickness of the growth face (Figure ge) . When the 
thickness is less than I** , the asymmetry of the probabil- 
ity distribution is reversed (Figure ||f). It is, therefore, 
through the combined action of the two thermodynamic 



constraints on the stem length that the thickness con- 
verges to a value close to /min- 




FIG. 7. Cut through a polymer crystal which was pro- 
duced by the growth of twenty successive layers on a growth 
face with a uniform thickness of 50 units at T = 2.0efc~^. 
The stems are represented by vertical cuboids. The cut is 16 
stems wide. 
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FIG. 8. Probability distributions of the stem length for a new crystalline layer grown at T = 2.75 ek~^ on a growth face 
of thickness: (a) oo, (b) 100, (c) 50, (d) 25, (e) 21 and (f) 19. The dashed vertical lines in the probability distributions are 
at the thickness of the growth face. 



The picture is not quite this simple at ah tempera- 
tures. As the supercoohng decreases, it becomes increas- 
ingly unfavourable for a stem to overhang the edge of 
the growth face. Indeed, for sufficiently small supercool- 
ings the probability distribution for the stem length never 
becomes symmetrical about the thickness of the growth 
face, not even when the thickness of the growth face is 
close to /mill- This situation is illustrated in Figure |6|b. 
After the growth of two layers on a 50-unit thick surface, 
the crystal stops growing because the outer layer is too 
thin for a new layer to form. For these supercoolings, 
as in the SG model, the rounding of the crystal profile 
inhibits growth. To overcome this barrier requires a co- 
operative mechanism whereby a new layer takes advan- 
tage of (and then locks in) dynamic fluctuations in the 
outer layer to larger thickness. However, unlike the SG 
model, the current model has no interlayer dynamics — 
we attempt to grow a new layer on an outer layer that is 
static — and so growth stops. Despite this it is clear that 
if this interlayer dynamics could be included, it would 
again lead to steady-state growth close to l^i^i- 
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FIG. 9. A schematic picture of a two-dimensional slice 
(perpendicular to the growth face) through a lamellar poly- 
mer crystal which forms the basis of the two-dimensional ver- 
sion of the Sadler-Gilmer model. The three possible changes 
in configuration allowed by the model are shown (the dashed 
lines represent the outline of the possible new configura- 
tions). 

We should note that this cessation of growth was 
also found in_.the model of Frank and Tosi at low 
supercoolings.E3 Lauritzen and Passaglia were also aware 
of this effect, but they introduced an ad-hoc energetic 
term in their rate constants to prevent it.o However, in 
the restricted equilibrium model of Price this effect was 
absent.Eil In this study each new layer, but not the crys- 



tal as a whole, was allowed to reach equilibrium and so 
the kinetic constraint on the stem length is absent. 

Finally, we should note that our multi-pathway model 
is not parameter-free, and that, like most other mod- 
els of polymer GT'y.stallization (including the LH theorya, 
the SGL pipdeSEZlH and the earlier multi-pathway 
modelscj^EllLS) , for some choices of parameters (not those 
used here) the lamellar thickness begins to increase at 
sufficiently large supercoolingE3 This effect occurs be- 
cause the large driving force for crystallization at large 
supercoolings reduces the effect that the thickness of the 
growth face has in constraining the stem lengths. 



IV. THE SADLER-GILMER MODEL 

In this section we re-examine the model used by Sadler 
and Gilmer in order to see whether the mechanism of 
thickness selection that we found in the previous sec- 
tion for our multi-pathway model also occurs in the SG 
model. Sadler and Gilmer interpreted this model in terms 
of an entropic barrier. In particular, they argued that the 
rounding of the crystal profile gives rise to an entropic 
barrier, which can only be surmounted by a fluctuation 
to a squarer profile before growth can continue. As this 
barrier increases with lamellar thickness it constrains the 
thickness to a value close to ^min- However, we shall not 
dwell on this interpretation here, but instead direct the 
interested reader to a critique of this argument in Ref. 
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FIG. 10. The dependence of the thickness of a layer in 
the bulk of the crystal on the thickness of the previous layer 
at T = O.OSTm. The dotted arrowed lines show the thickness 
converging to the fixed point of the attractor from above and 
below. 

In the SG model the growth of a polymer crystal re- 
sults from the attachment and detachment of polymer 
units at the growth face. The rules that govern the sites 
at which these processes can occur are designed to mimic 
the effects of the chain connectivity. In the original three- 
dimensional version of the model, under many conditions 



the growth face is rough and the correlations between 
stems-JB^the direction parallel to the growth face are 
weakJlj'Ea Therefore, an even simpler two-dimensional 
version of the model was developed in which lateral cor- 
relations are neglected entirely, and only a slice through 
the polym8&,(pjystal perpendicular to the growth face is 
considered O'Ej 
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FIG. 11. (a) Probability distributions of the stem length 
in the bulk of the crystal given that the previous layer has 
the labelled thickness, (b) A contour plot of the log of a se- 
ries of such probability distributions. (Only thirty contours 
are displayed. The probability continues to fall rapidly in 
blank top left-hand corner.) The fixed-point attractor has 
been overlaid on this figure to illustrate its connection to the 
probability distributions. T — 0.95Tm- 

The geometry of the model is shown in Figure M. 
Changes in configuration can only occur at the outermost 
stem and stems behind the growth face are 'pinned' be- 
cause of the chain connectivity. At each step, there are 
three possible changes in configuration: the outermost 
stem can increase in length, a new stem can be initiated 
and a polymer unit can be removed from the outermost 
stem. The model can be formulated in terms of a set 



of rate equations that can be easily solved by numerical 
integrationJid 

When we examine the dependence of the thickness of 
a layer on the previous, we again find a fixed-point at- 
tractor describing the convergence of the thickness to its 
steady-state value (Figure |l^). Moreover, when we ex- 
amine the probability distributions for the stem length 
we find evidence for the same three constraints as for the 
multi-pathway model (Figure |llp). The weaker nature of 
the kinetic constraint is particularly clear from the much 
more rapid exponential decay of the probability for stems 
that extend beyond the growth face. The role played by 
the two thermodynamic constraints in the mechanism of 
thickness selection is particularly clear from Figure Ob. 
As the thickness of the growth face decreases the viable 
range of stem lengths decreases until the the thickness of 
the growth face meets Zmin at the fixed point. 



V. DISCUSSION 

In this paper we have outlined evidence from com- 
puter simulations for a mechanism of thickness selection 
in lamellar polymer crystal that differs from the theories 
of Lauritzen and Hoffman, and Sadler and Gilmer. In- 
stead, the mechanism has much more in comjaon .with 
the results of earlier multi-pathway modelsE3'ElrE3 We 
find that a fixed-point attractor which describes the dy- 
namical convergence of the crystal thickness to a value 
just larger than the minimum stable thickness, ^min- This 
convergence arises from the combined effect of two con- 
straints on the length of stems in a layer: it is un- 
favourable for a stem to be shorter than /min and for 
a stem to overhang the edge of the previous layer. It is 
encouraging to note that we find the same mechanism of 
thickness selection operating in two models which make 
very different assunrptions about the microscopic growth 
processes. This provides evidence of the generality of this 
mechanism, and so suggests that, although the models 
described here have a very simplified description of the 
microscopic dynamics, the physical principles behind the 
mechanism could be general enough to apply to real poly- 
mers. 

This mechanism of thickness selection is also consis- 
tent with experimentS|-wiiere the temperature is changed 
during crystallization. c3^EZl The steps that result indicate 
that the thickness of the lamellar crystals dynamically 
converges to the steady-state thickness for the new tem- 
perature by a mechanism similar to that which we ob- 
serve in our simulations. Furthermore, if the step pro- 
files could be characterized with sufficient resolution by 
atomic-force microscopy, it may be possible to extract the 
fixed-point attractor of a real polymers. However, for a 
temperature decrease the step profiles may also refiect 
the rounding of the crystal edge and foi^ temperature 
increase the roughness of the fold surface.cJ Furthermore, 
any annealing mechanisms that operate could change the 



shape of the step profile from its as-formed state. 

Although the multi-pathway approach is, in some 
ways, an extension of the LH theory, the removal of many 
of the LH constraints leads to significantly different be- 
haviour. In particular, our work undermines the LH as- 
sumptions that the initial nucleus determines the thick- 
ness of a layer, and shows that the approach embodied in 
Equation (g) (i.e. a comparison of the growth rates of the 
crystals in an ensemble of crystals of different thickness 
all of which grow at constant thickness) is inappropriate 
because crystals of arbitrary thickness do not necessarily 
continue to grow at that thickness. Although our results 
lead us to question the thickness selection mechanism in 
the LH theory, other aspects of the nucleation approach 
may not be affected by our critique. For example, the 
regime transitions are a result of the different functional 
dependence of the growth rate on the nucleation rate a: 
the substrate completion rate in the different regimes. 

Recently, there have been a number of alternative the- 
oretical proposals that have made recourse to metastable 
phases. Keller and coworkers suggested that crystalliza- 
tion of polyethylene could initially occur into the mo- 
bile hexagonal phase. These crystals would then thicken 
until a critical thickness was reached at which a ptf^ 
transition to the orthorhombic phase would occur Ji3^EJ 
Olmsted et al. have argued that the density fluctua- 
tions resulting from the spinodal decomposition of a poly- 
mer meltc3 assist the nucleation of crystals.c3 Strobl and 
coworkers have argued, on the basis of the thickness de- 
pendence of the crystallization and melting temperatures 
of syndiotactic polypropylene, and the granular texture 
in AFM images of the same polymer, that the polymer 
first crystallizes into blocks, which ar|a^ubscqucntly sta- 
bilized when they fuse into lamellaeo Our simulations 
can say little about these proposals since our polymer 
models are too simple to be able to capture such fea- 
tures. However, all these approaches are based on be- 
haviour that has been observed in crystallization from 
the melt, so it is not clear how the ideas can apply to 
crystallization from solution, where the same basic laws 
for lamellar polymer crystals apply. 
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